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Abstract 

We study the spectra of scalar and vector mesons in four dimensional strongly coupled SQCD-like 
theories in the Veneziano limit. The gauge theories describe the low energy dynamics of intersecting 
D3 and D7-branes on the singular and deformed conifold and their strong coupling regime can be 
explored by means of dual fully backreacted supergravity backgrounds. The mesons we focus on are 
dual to fluctuations of the worldvolume gauge field on a probe D7-brane in these backgrounds. As 
we will comment in detail, the general occurrence of various UV pathologies in the D3-D7 set-ups 
under study, forces us to adapt the standard holographic recipes to theories with intrinsic cutoffs. 
Just as for QED, the low energy spectra for mesonic-like bound states will be consistent and largely 
independent of the UV cutoffs. We will study in detail how these spectra vary with the number of 
the fundamental sea flavors and their mass. 
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1 Introduction 



Holography is nowadays a standard and powerful method to investigate properties of some 
strongly coupled gauge theories. Most of the studies of flavor physics in this context have 
focused on the quenched approximation where the internal quark loops are neglected. The 
extension of these investigations to the unquenched cases, where the full dynamics of the 
flavors is included, is of obvious interest. The main focus of this paper is the study of mass 
spectra of low spin (J = 0, 1) mesons in certain strongly coupled SQCD-like theories in the 
unquenched Veneziano regime, where the number of colors and the number of flavors are 
both taken to be very large, with their ratio taken to be fixed. The aim is to extract the 
dependence of the spectra on the number, Nf, of sea flavors and on their mass, rriq. 

The Af = 1 SQCD-like models we will consider describe the low energy dynamics at the 
4d intersection of "color" D3-branes and Nf homogeneously smeared "flavor" D7-branes 
on the singular and on the deformed conifold. The models are the flavored unquenched 
versions of the conformal Klebanov-Witten (KW) [Ij and the confining Klebanov-Strassler 
(KS) [2j ones. The strong 't Hooft coupling regime of these theories is mapped to some dual 
backgrounds, arising as the full backreaction of the color and the flavor D-branes, with Nf 
D7-brane sources. The supergravity background dual to the KW (resp. KS) model coupled 
to chiral (resp. non chiral) massless dynamical flavors was found in [3] (resp. [1]). The 
solution with chiral massless flavors in the KS case, which we will not consider, appears 
in [5]. These massless-flavored solutions have generically (good) singularities at the origin 
of the transverse radial coordinate. The singularity is avoided in the massive case. The 
supergravity dual of the KW model coupled to chiral (resp. non chiral) dynamical massive 
flavors was found in [6] (resp. [7]). The relevant solution in the KS case with non chiral 
massive flavors was found in [8]. 

In order to extract the mesonic spectra we will probe these backgrounds with an external 
D7-brane and study the fluctuations of a selected set of decoupled modes (dual to scalar 
and vector mesons) of the gauge field on its worldvolume. This study will require a careful 
treatment of the boundary conditions of the fluctuating fields, in connection with the fact that 
the backreacted D3-D7 backgrounds - unlike the corresponding unflavored ones - have various 
UV pathologies. In particular, we will show that there is a singularity in the holographic a- 
function at finite radial position. Just as in QED, which has a UV Landau pole, the spectra 
of the various bound states will nevertheless be meaningful and independent of the cutoffs 
(up to corrections suppressed by the UV scale). 

Since the goal of this paper is to study the dependence of the mesonic spectra on the 
number of sea flavors and their masses, we will have to decide how to interpret our results, 
i.e. how to compare different theories with different flavor parameters. As already noticed 
in [8], there is no obvious natural scale or coupling which can be assumed to stay fixed when 
varying those parameters. In any case, the comparison between two theories will necessarily 
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require us to make a choice of what to keep fixed. The comparison will strongly depend 
on this choice. To disentangle the possible different conclusions we will need some explicit 
formula for the mesonic masses as a function of the physical parameters in the theory. 

We will parameterize the meson masses as Mmeson = ^ where Aq is a dimensionful scale 
which can be "measured" in the IR and which we will take fixed when comparing different 
theories. This scale will be related to the "coupling" and the effective "string tension" at 
the probe quark mass scale. The coefficients uj will then be determined from a numerical 
analysis. 

Our results show that the tu's receive small "radiative" corrections due to dynamical fiavor 
loops and decrease as the effective coupling gprpNj (where g^j. ~ e*^) is increased (or if the 
mass of the sea flavors is decreased). These small "quantum effect" corrections resemble the 
hydrogen atom Lamb shift in QED. In section [5] we will argue that, for the flavored KW 
models, these corrections could be related to the small variation of the effective number of 
(adjoint plus bifundamental) degrees of freedom due to the internal flavor loops. In the KS 
cases they could be related to corrections to the glueball (or Kaluza Klein) mass scale. The 
decrease of the cj's depends on the choice of keeping flxed the coupling at the probe quark 
mass: had we kept flxed the coupling at some larger UV scale, larger gprpNf could have 
yielded larger meson masses. 

The structure of the paper is the following. We will start, in section [21 by examining the 
general UV behavior of the D3-D7 models under study. We will show how non trivial UV 
cutoffs emerge when the holographic a-function is considered. Moreover, we will outline the 
general limits of our analysis. We will then focus on scalar and vector mesons in the flavored 
conifold models. In section [3], combining the numerical shooting technique and the WKB 
analysis, we will extract the spectra of mesons having massive constituent flavor flelds in 
the (flavored) KW models. Having chosen a prescription to compare different theories we 
will study how the sea flavors (either massless or massive) affect the mesonic spectra. In 
section H] we will present the results of analogous studies in the flavored KS models, limiting 
our analysis to the case of mesons with massless fundamental constituents. In section [5] we 
will discuss a possible way to interpret our results. We will end in section [6] with a set of 
concluding remarks. Various technical details and some comments on high spin mesons will 
be left to the appendices. 

2 General comments 

A common feature of D3-D7 set-ups, both in the case of smeared and localized [9] D7-branes, 
is a running dilaton in the dual supergravity solutions. The dilaton, and consequently the 
effective string coupling, typically increases with the transverse radial variable u of the 
backgrounds and blows up at a certain point Ud- The solutions can thus be deflned at best 
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only in a region u < Ucut < Ud where the dilaton can stay smallo 

The point Ud can be formally mapped to a UV Landau-pole in the dual field theories. 
The simplest way to realize it is to consider a localized D3-D7 set-up in fiat space. The 
theory on N^. D3-branes is just A/" = 4 SYM, which has exactly zero beta function. We can 
add Nf fiavors to the theory by means of a stack of D7-branes, sharing with the D3-branes 
the 4d Minkowski directions and extended along a non compact 4-dimensional submanifold 
of the transverse space. This breaks supersymmetry (to A/" = 2) as well as conformal 
invariance. The perturbatively exact coefficient of the beta function for the inverse squared 
gauge coupling is 60 = SAi'c — 'iNc — Nf = —Nf, which means that the theory has a UV 
Landau pole. This pathological UV behavior is inherited by D3-D7 set-ups both in orbifold 
and conifold models. 

Apart from the dilaton divergence, the theories we are going to consider have other less 
explicit UV singularities. In the fiavored KW cases the integration constants can be fixed so 
that the warp factor vanishes exactly at the Landau pole Ud- In the fiavored KS cases this 
is not a consistent choice and the warp factor will diverge badly (going to minus infinity) at 
Ud- This means that there will be a point Uh < Ud where the warp factor vanishes and the 
metric becomes singular. A possible interpretation of this point is in terms of a duality wall 
in the dual theory |1] . Since for u> Uh the warp factor becomes negative, Uh replaces Ud as 
a sensible UV "end" of the fiavored KS backgrounds. 

Both the fiavored KW and KS models have an even more "hidden" potential UV pathology. 
It shows up by considering their 5d reduction. Starting with a lOd string frame metric of 
the form 

ds^ = a{u) \dx^dx^ + I3{u)du^]^ + dsl^^ , (2.1) 

the reduction on the internal five dimensional manifold (with volume Vint{u)), gives a 5d 
Einstein frame metric of the form 

dsl = H{uf/^ [dx^dx^ + (3{u)du'^] , (2.2) 

where 

In standard set-ups, the function H(uY^^, which can be roughly identified with the dual field 
theory energy scale, monotonically varies with the radial coordinate. This is also required 
in order for the "holographic a-function" |10] 

a{u) ~ /5(u)3/2if(M)^/2 [//'(«)] -3 ^ (2.4) 

^ Just to fix the notation consistently with the following sections, we will call p (resp t) the radial variable 
in the flavored KW (resp. KS) solutions and we will set Ud = Plp (resp. Ud = Td). In our conventions large 
values of u are mapped to large energy scales in the dual field theory. 
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Figure 1: The function H in the massless-fiavored KW modeL 

to be finiteil 

In both the flavored KW and KS cases, instead, the function H{u) is not monotonic: it 
increases with u from zero up to a maximum at a point Ua and then it decreases back to zero 
(at Ud in the KW cases and at Uh, with Ua < Uh < Ud in the KS cases). A representative plot 
is given in figure [H This behavior implies that the holographic a-function a{u) as defined in 
(12. 4p is singular and discontinuous at Ua- As we will show in the following, the a-function 
will be positive and increasing up to m = u~ where it will asymptotically diverge. At u = 
a{u) will go to minus infinity and then it will increase to zero (at Ud in the KW case and 
at Uh in the KS one) for increasing u. In order to avoid the apparently pathological region 
Ua <: u < Ud, we will cutoff our D3-D7 backgrounds at Ucut < Ua- 

The occurrence of UV pathologies in the D3-D7 models we are going to consider, does 
not subtract interest to the study of their IR dynamics. QED, the most notable quantum 
field theory model with a UV Landau pole, certainly provides excellent and meaningful 
predictions on the low energy (e.g. atomic) spectra, which do not sensitively depend on the 
UV cutoff. In order to study, say, the spectrum of bound states of (s) quarks pairs in our 
D3-D7 conifold models we will just have to reconsider the standard holographic recipes and 
adapt them to the case where the dual supergravity backgrounds do not have a boundary 
at infinity. Having discussed the presence of UV pathologies, we stress that the aim of this 
paper is not to study how to cure them or how to consistently UV complete the theory. 
Conversely, we focus in computing IR quantities that are only mildly affected by whatever 
the UV physics is, in a sense that we will make precise below. 

In the following we will focus on vector and scalar mesons dual to fluctuations of the gauge 
field on the worldvolume of a D7-brane probe. The corresponding mass spectrum will be 
quantized after imposing that the fluctuations be regular at uq (the minimal radial distance 
reached by the brane) and vanish at the cutoff Ucut < Ua- The pathological UV region 
will thus be excluded and treated as producing an infinite wall in the effective quantum 

^Thc monotonicity of H{u) also plays a crucial role in holographic computations of the entanglement 
entropy, see jllj . The notations of that paper are used in the equations above. 
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mechanical description of the fluctuating modes. The choice of the cutoff (but not of its 
maximal value) is nevertheless arbitrary in the range uq <^ Ucut < Ua and the consistency of 
our results will be guaranteed provided we show that the dependence on the cutoff is highly 
suppressed. Note that we can think about our models as the IR regions of UV consistent 
theories, providing some UV completions to the backgrounds at hand. This perspective has 
been employed recently in ^12i| . 

Let us conclude with the following remarks. The theories to which we add flavors cor- 
respond to the low energy dynamics of regular (and M fractional) D3-branes on the 
(deformed) conifold. They are (cascading) A/" = 1 4d gauge theories with gauge group 
SU{N) X SU{N) {SU{N + M) X SU{N)) and bifundamental matter flelds A,B trans- 
forming as SU{2) X SU{2) doublets and interacting with a quartic superpotential Wkw = 
e^H^^ AiBkAjBi. The perturbative superpotential in the flavored case is: 

W = Wkw + hi gi(Ai5i - Aa^a)?! + /12 q2iBiAi - 53^2)^2 + h {mif + m {mi) , (2-5) 

where qi and are the flavor multiplets and hi and ki are the couplings. 

The SQCD-like models under study are assumed to be in the Veneziano regime. This 
means that we take Nc,Nf 00 (where Nc is the number of colors) with Nf/Nc and 
A = g^prpNc flxed. Moreover, the dual supergravity solutions (with DBI-I-WZ source terms 
for the smeared D7-branes) will only be reliable provided A ^ 1 (as usual) and Nf <ti 
These limits will allow us to consider cases where the effective coupling g%rpNf = NfX/Nc, 
weighting the vacuum polarization effects due to the dynamical flavors, is of order one|j 



3 Meson excitations in the KW models with flavor 

In this section we will discuss (part of) the spectrum of meson masses in the particular 
framework of the so-called Klebanov-Witten conformal model [1] and its generalizations with 
massless |3j and massive |7j unquenched flavors. We start by writing general expressions for 
the background and excitations and then specialize the study for the different cases. We will 
close the discussion of the flavored KW models by analyzing the holographic a-function in 
section 13.61 

■^It is worth pointing out that the Nf <^ Nc condition, necessary in all the D3-D7 models, is not mandatory 
in different set-ups with other kinds of brane intersections. Some examples where Nf/Nc can be kept of 
order one have been studied in [13l [TH [15] . 

"'Note that the open string coupling on the flavor branes is not Nfc'^: the backreacting branes are smeared 
in the transverse space, so only a small fraction of them is within a distance of order \/a' and the coupling 
is parametrically smaller than Nfc"^ [51 115j. 
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3.1 The background solution and the excitation equations 



We consider solutions of type IIB supergravity coupled to a homogeneously smeared set of 
Nf D7-branes (we refer the reader to [31 El Ej for further details). The ansatz for the fields 
which take non-trivial values in the solution is, in Einstein frame: 



dsl 



^(5) 



^(1) 



h 2dxl;^ + a'h2 [e^^ dp^ + dsl] ' 
— 2^ {de^ + sin^ Oidifi'^i) + ^{di^ + cos Oidip 



i=l,2 



2 



i=l,2 



d'^x A d{h ^) - '^^^-^"'^ sin sin Q^dOi A difi A d92 A d(p2 A dijj 



9sNf{p) 



{dip + cos 9id(pi + cos 62dip2) 



Notice that the 5-form is self-dual and, using the relations: 

11 „ 1 



2k 



T 



gs{2Ti)Pa'— 



(3.i: 



(3.2) 



one can check that it satisfies the quantization condition J^__ F(5) = A^c^sSk^^^q-j = (27r)^a'^(7sA''c, 
where is the internal manifold. If the two-form VL is the density distribution of the smeared 
D7-branes {Qs^ = — c?-F(i)), the action of the gravity+branes system reads ^: 



S 



R - -Omc 



^(5)1 



+ 



(3.3) 



Then (13. ip provides a supersymmetric solution if: 



3^?siV/(p) ^^ 
47r 



/ = 3 - 2e2/ 

/i = —27'KgsNce 



2g _ ^9sNf{p) ^^ 

8tt 



(3.4) 



These equations are valid for the case without unquenched flavors [T] {Nf{p) = such that 
g = f = p), the case with massless unquenched flavors [3J (A^/(p) = const > 0) and the 
case with massive unquenched flavors |3l[6l[7] (where iV/(p) becomes a suitable p-dependent 
expression depending on the kind of flavor branes under consideration). 

In the following, we analyze excitations of a brane probe which preserves the same su- 
persymmetry as the background. Thus, there are two kinds of flavor branes in the generic 
set-up: the first kind corresponds to dynamical quarks, i.e. the Nf branes backreacting on 
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the geometry, accordingly corresponding to unquenched flavors (in section [3^ we will denote 
Pq the value of p at their tip). The second kind is a single oscillating probe brane (with tip at 
P — Pq)i associated to the quarks which actually constitute the mesons under consideration. 
Of course it is possible to take this single fluctuating brane to be one out of the Nf back- 
reacting ones but in the following we deal with the more generic case of allowing non-equal 
masses for dynamical and test quarks (p^ 7^ pq). Since the flavor branes affect the values of 
the metric, dilaton and p-forms in the solutions, one may wonder whether it makes sense to 
consider the oscillation of a brane in a background where the closed string fields are taken 
to be constant. That is indeed the case: a meson is associated to the oscillation of one (or 
a pair of) flavor brane and not to a collective oscillation of a set of order Nf flavor branes. 
Thus, the possible backreaction of the mesonic oscillation on the closed string background 
strictly vanishes in the Veneziano limit. 

Let us consider a flavor brane probe of the type first discussed in |16j and corresponding 
to the embedding 2:4 — ^3 = e2 where ^3,4 are two of the complex coordinates defining 
the conifold. Details and notations are spelled out in appendix lA.li We will only study a 
reduced subset of all possible mesonic modes. Concretely, we will not discuss fluctuations 
of the embedding and only the following fluctuations of the worldvolume gauge field, which 
correspond to a vector and a scalar in the dual gauge theory: 



Here, is a constant transverse vector and hi is the (angular) left-invariant one-form defined 
in flA.91) . We have chosen this particular set of fluctuations for the sake of simplicity, since, 
as we will see below, it gives rise to relatively simple decoupled differential equations both 
in the present set-up and in the Klebanov-Strassler case to be discussed in section HI Also 
for simplicity, we have not included any angular dependence for the fluctuating gauge fleld 
(a truncation which is non-trivially consistent). Even if we restrict ourselves to this very 
limited subset of all the possible excitations, we expect that the results can point out the 
general trends on how the presence of dynamical flavors affects the meson masses. 

The procedure to obtain the second order equations associated to these modes is outlined 
in appendix IA.2[ One gets: 



A = a^{p)e'^''i^dx^' + a,(p)e*^^/ii . 



(3.5) 







dp {e^s-^'Pie^P - e'P^)dpa:) + Mla'he^^+^f (l + e^^^-^'Pi-^e^^-^f - 1) j a, , (3.6) 







'o I ; dgas I dnih ^) — 5 5-— 





(3.7) 



^An alternative embedding has been considered in [T7] . 
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where M^, = — A;^ are the masses for the vector and scalar mesonic excitations under con- 
sideration. As it is customary, enforcing appropriate IR and UV behaviors for the functions 
av{p),as{p) will select a discrete spectrum for My^s- 

Before turning to the study of this issue in several cases, let us introduce a useful param- 
eterization of the meson masses, factoring out a dimensionful scale. Define: 



(3.8) 



where Tq is the tension of a hyp_othetical 
and Xq is the 't Hooft coupling 



fundamental string stretched at constant p 
^ at the same scale p = pq: 



PQ 



Tr 



>/2 



a/ GttGxx 



\P=PQ ' 



A. 



87ig,N,e^ 



\P=PQ 



(3.9) 



The important point is that these are quantities "measured" in the IR, i. e. at the tip of the 
brane p = pq, which is related to the quark mass. The u^^s "will be towers of numbers which, 
as we will see, depend on j^Xq. 

Notice, however, that, when writing (13. 8p . we could have chosen to factor out a different 
IR dimensionful scale. Basically, this is related to the observation in [18] of the possibility of 
different definitions of the radius-energy relation in non-conformal theories. For instance, we 
could have factored out the constituent quark mass to be discussed in section 13. 5[ Thus, one 

should keep in mind that when we compute how u varies with jfXq, we will be computing 

1 1. 

how meson masses change while keeping Tq / Xq fixed. 



3.2 Mesons in the quenched approximation 

We want to analyze the equations (13. 6p . (13.71) when the background is just the unflavored 
KW solution, i.e.: 

27 

Nf{p) = 0, f = g = p, = const , h = —TxgsN,e-^P . (3.10) 

It is useful to define: 

P = P-PQ- (3-11) 

^ The value of the tension Tq of eq. (13. 9|) can be obtamed by studying the short distance behavior of 
the QQ potential, obtained by analyzing a hanging open string in the background metric, see appendix IdI 
Notice that Tq is simply the fundamental string tension l/(27ra') redshifted by the string frame warp factor. 
For the identification with the 't Hooft coupling, we have used the orbifold relation ^ + = g~^e~'^ . 

Let us consider g\ — g\ = g\rj^ = Sngge'^, and thus define the 't Hooft coupling as A = Ncgprp — SirggNce'^. 
One should keep in mind that this orbifold relation is not guaranteed to hold in general. The relations (|3.9p 
should be taken as the definition of Tq, Xq. 
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The equations (13.61) . (13.71) then read: 



e-3p 



= dp{e-\e^~'-l)d-Mp))+^l{^-^Xp) . (3.12) 

_33(=3^ — 94 _ p~'^P 

= d,{e\e'-p - l)dMP)) - _ ^) ^s{fi) + ^y~'i\ + ^)a.(p) , (3.13) 

where we have used the definitions (13.81) and (13. 9p . In this conformal case, the actual quark 
bare mass coincides with the constituent mass defined as the energy of a string stretched 
from the bottom of the geometry up to the brane tip pg: 

1 rPQ ^ (>%(>PQ 
ran = / Va'e^e^dp = = . (3.14) 

This can be used to rephrase (13.81) as: 

M,,, = J^u;,,, . (3.15) 

Comparing (13. 8p to (13.151) , we see that Tq ~ mg / a/Aq . 

We can now numerically analyze equations (I3.12p . (I3.13p . In each case, requiring regularity 
at p = and normalizability in the UV selects a discrete set of values for the cu's. A standard 
computation using the shooting technique yields: 

LUy = 2.337, 4.720, 7.088, 9.454, . . . 

oj, = 5.174, 7.358, 9.482, 11.580, ... (3.16) 

3.2.1 The Schrodinger potential formalism and WKB estimates 

As shown in appendix [B], after performing a convenient change of variables, the fluctuations 
equations (I3.12p and (I3.13P can be written as Schrodinger equations for some particular 
potentials. We can then apply the WKB approximation to get an estimate of the mass 
levels. In the case of the fluctuations described by eqs. (I3.12p and (13.130 this analysis leads 
to the following estimates of Uy^. 



/ 7 15 

4")^C.^, + + - , (n = l,2,...) (3.17) 
where and are given by: 

Co ~ 2.365 , Cs ~ 2.051 . (3.18) 

In fact, from WKB we only get terms of order and n° as n — > oo. Thus, the term y 
inside the square root of the expression for does not come from WKB and we obtained 
it by fltting the numerical data. In flgure [21 we compare these expressions to the flrst few 
numerically found eigenvalues. 
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Figure 2: A plot of the first few values of Uv and Us obtained numerically (dots) and the 
approximation given in eq (13.171) (solid line). 

3.3 Mesons in the presence of unquenched massless flavors 

In this section, the background defined by the equations in (13.41) with constant Nf{p) = 
Nf > [5] will be considered. The relevant solution is: 

e = 

3gsNf{pLp - p) 

= Cs^/QipLP - P) (1 + QiPLP - p))-^e(^-^^^) , 
= c3(l + 6(pLP-p))^e(^-^^^), 

h ^ ™i,(J^)^(r(i,-^-4(p..-p))-r(i,^ , (3.19) 

where the F represents the incomplete gamma-function. 

It is worth commenting on the integration constants appearing in (I3.19p . As compared 
to [3j, we have explicitly kept the integration constant pip which fixes the position of the 
Landau pole (where the dilaton diverges). We have also kept the constant C3 which can be 
reabsorbed by rescaling the Minkowski coordinates and thus just rescales what one defines as 
energy. On the other hand, we have set the constant ci (see eqs. (2.37), (2.38) of [3j) to zero 
for the sake of IR regularity. Even if the unquenched solution with massless flavors is always 
IR singular, it was shown in [3] that ci = produces the less severe IR singularity. Most 
importantly, it was explicitly shown in [Gj [7] that by introducing any non-zero mass for the 
dynamical quarks, a regular IR can be obtained. Taking the massless limit of this family of 
massive regular solutions yields the ci = condition. Finally, we have set /i(pLp) = 0, such 
that metric and dilaton are both singular at the same point p = pip. This UV prescription 
is not important for the computation of the meson masses: enforcing, instead, h{ph) = 
for some ph < Plp would just add an additive constant to h. As long as ph ^ Pg, this 
would only modify the values of the meson masses by quantities exponentially suppressed as 
gPQ-Ph Extended comments regarding the effects of UV prescriptions on the computation 
of the masses can be found in sections 13.3.11 and 13.3.21 
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By inserting (13.191) into (13. 6p . (13. 7p . and using again the definitions (13.81) . (13.91) . one finds 
the appropriate second order equations. It is easy to check that the resuhing spectrum of 
id's only depends on p^p — pg, which can be related to the physical quantity jfXq by using 
the explicit solution for the dilaton in (13. 19^ : 



The quenched case is recovered when j^Xq — or Plp — Pq — ^ oo. 
3.3.1 Schrodinger potentials and UV cutoffs 

Our goal now is to study equations (13. 6p . (13.70 in this background. In each equation we have 
to demand regularity at pq for the function that defines the fluctuation. On the other hand, 
one cannot use the usual UV normalizability condition since there is a Landau pole which 
severely modifies the UV behavior. As anticipated in section |2l the natural condition is to 
require for each excitation that: 



These conditions amount to introducing by hand an infinite wall at pcut- Our working 
prescription will be to set pcut ^ Pa, i-^- we consider a UV completion of the theory slightly 
below the scale in which the holographic a-function becomes singular. 

We now turn to justifying the proviso (I3.2ip . Schrodinger potentials associated with 
the fluctuation equations can be computed following the definitions of appendix [Bl Figure 
|3] shows two examples of Schrodinger potentials for the vector excitation. For the scalar 
excitation, the plots are qualitatively similar (except that V diverges towards +oo in the 
UV pathological region). 

From the figure, it is apparent that as long as Plp is large {jfXq is small), the UV 
pathological region lies very far from the minimum of the potential where the physical wave 
function has its main support (see the plot on the left). On general grounds, we expect 
that whatever would be the UV completion of the potential, it only affects the eigenvalues 
by exponentially suppressed quantities. In particular, the precise value of pcut < pa < Plp 
has only a negligible effect on the spectrum, as long as pcut lies far away from the IR region 
where the potential reaches its minimum. 

On the other hand, if j^Xq is too large, pip (and thus Pa, Pcut) is not far from the IR 
region. Then, the formalism breaks down. In physical terms, if the UV completion sets in 
not far from the relevant IR scale, it will affect the IR physics in a non-negligible way. The 
plot on the right of figure [3] shows an example of this behavior. For this reason, we will 
restrict ourselves to jfXg < 10. 



Plp = Plp - Pq 



SNfXQ 



(3.20) 



aviPcut) = asipcut) = . 



(3.21) 
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Figure 3: On the left, the Schrodinger potential for the vector excitation with = 3, 
pLP = 20. The UV pathological region is well separated from the IR region. On the right, 
the potential for = 3, plp = 2. Obviously, the UV region is not separated from the IR 
one. The potentials diverge to minus infinity ai y = log{e^^^ — 1). 



3.3.2 Estimates of the meson spectrum from WKB 

We have checked that the WKB expressions (13.171) in the unquenched case are also in very 
good agreement with the numerical values obtained from the shooting technique|l| Thus, we 
just concentrate on studying the WKB integrals (lB.16p . (IB. 170 . 

For small NfXg/Nc (large pip — Pq), we can expand the integrands in flB.16|) . flB.17|) . A 
simple computation shows that: 

C ~ 2.365(1 - 1.12 X 10-=^-^ Aq + ...) , 

C ^ 2.051(1 - 1.08 X 10"^-^ Aq + ...) . (3.22) 

It is clear from (13.221) that the leading correction to the masses, due to the dynamical flavors, 
is small. It is also approximately linear in NfXq/Nc = Nfgpj,\p=pg, see the bottom lines in 
figure m The fact that the meson masses decrease with Nf might appear counter-intuitive 
at first sight. In fact, there is nothing strange about it, as we will argue in the subsection 
13.3.31 Moreover, a qualitatively similar behavior can be inferred from the unquenched lattice 
calculation reported in [19] . 

Now that we have found (13.221) . let us be more precise about the size of the uncertainties 
of the spectrum associated to UV prescriptions. For large p (but p < Plp), the integrand of 
the WKB integrals (IB. 161) . (IB. 171) behaves as e~^. Thus, if one takes a different value of pcut, 
the variation of the integral and thus of the masses is of order e~^'="*. Taking into account 
the radius energy relation e'' ~ A (where A is the energy scale), one sees that different UV 



^The relative error of the (|3.17p formula with respect to the obtained numerical values is always well 
below 1%, except for uji'^\ when in the quenched case the error is already 1.2% and in the unquenched cases 
remains approximately of the same order. The same applies to the case of unquenched massive flavors to be 
studied in the next section. 
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completions yield variations of the meson masses of order ^f^, where the IR scale is here 
associated to the mass of the quark constituents of the meson and Auv the scale at which 
the UV completion sets in. Notice that these corrections are far more suppressed than the 
first corrections displayed in fl3.22p . which are of order {pLp)~^ (logarithmic in the energy 
scale) . 



3.3.3 Comments on the physical interpretation of the results 

From eqs. (13.81) . (I3.17p . (I3.22p . we read the following expression for the tower of vector 
mesons: 

N 

Mi") ^ 6.19 n^(l - 1.12 X IQ-^-IAq + . . . ) . (3.23) 

\A Nc 

From this formula, it seems that meson masses decrease with Nf, what, in turn, seems 
to contradict the fact that dynamical quarks screen the color charges: the binding energy 
between quark and anti-quark should decrease with increasing Nf and, accordingly, meson 
masses should increase. The crucial point is that the flavor effects on the meson masses 
heavily depend on the scales we keep fixed while comparing different theories. This means 
that (I3.23P states that meson masses decrease with Nf if the IR quantities Tq and Xq are 
kept fixed when comparing theories with different Nf. Had we kept fixed the gauge coupling 
at some UV scale, larger Nf would have yielded smaller Xq, possibly resulting in larger 
meson masses. We now clarify this statement with an example. 

Suppose we want to rewrite (13.230 in terms of Tq and of the 't Hooft coupling A,,, at a 
scale > Pq- Due to the Landau pole, it does not make sense to take in the far UV, and 
in fact let us assume p^ -C Plp- Using (13.91) . (I3.19p . (13.201) . we find that: 

A« = A.(1-(P.-P«)|^ + ...) . (3.24) 

Equation (I3.23P can be rewritten as: 

Mi") ^6.19 n -4 (^1 + 1. 12 X 10-3^A*(-l + 2.12(p, -pq)) +...^ . (3.25) 

Thus, when we compare theories with different Nf keeping fixed Tq and the coupling at any 
scale p^> Pq + 0.47, masses do indeed increase with Nf. 

Analogous considerations apply to the scalar mesons and to the rest of the cases discussed 
later in the paper as well. 
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3.4 Mesons in the presence of unquenched massive flavors 



We can generalize the analysis of the previous section to the case in which the backreacting 
dynamical flavors are massive. We take these dynamical flavors to be of the same (non-chiral) 
type as the probe one, but with different mass nig (such that the tip of the backreacting 
branes lies at Pg). The background geometry is obtained by inserting in (13. 4p the expression 
for Nf{p) computed in [7j, namely Nf{p) = for p < pg and Nf{p) = Nf{l — e^''«~^'') for 
P ^ Pq- Let us define: 



ki = l + 6{pLP -p)+ 2e^^'(e~^^^^ - 2e~^P) + e 
k2 = 6{plp - p) + 2e^P''{e-^P'^'' - e 
fc, = 2(e3''^-=^^^^-l) + 6(pip-p,) 



3pN 



(3.26) 



so the functions determining the background can be written as: 



57r 



9sNfk2 



e9 = C3 kY^e^'P^^ 
rpLP 

/ kf'\~^P+^P^Pdp 
J o 



.1/2 



03 47^6"-^^^ 



1/3' 



(P > Pq) 



(3.27) 



and: 



Svr 



QsNfkq 



e3 = ef = C3 k]l^eP-P^P 



'^3 \Jpq ^ 



2/3|'„-4p _ -4p, 



(P < P,) . (3.28) 



With this input, we can easily analyze the WKB integrals (lB.16p . (IB.lTp . Now, they depend 
on two parameters, namely and the mass of the dynamical quarks through the quantity 
Pq~ P^- Results are plotted in figure HI At fixed Pg — Pq, the meson masses shghtly decrease 
^ ^ As -rfXo 0, the C's tend to the quenched values (13.181) . Moreover, the larger 



with i^A, 



the value of pg — pg, the larger are the meson masses at fixed ^Aq. For large Pg — pg, 
the C's tend to the quenched values too, i.e. the lines in figure H] become horizontal. This 
is expected since if the dynamical quarks are very massive, they become quenched, even if 
there are many of them. At small Pg — pq — oo, one recovers the massless quark set-up of 
section 13.31 and indeed the lower line in each plot is well approximated by (I3.22p . 



3.5 On the constituent quark mass 

As explained above, the definition (13.81) . which uses the local string tension, and (I3.15p . which 
involves the constituent mass defined as in (13.141) . coincide in the quenched (conformal) case. 
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Figure 4: The functions Ci),s entering the meson mass formulae, see eqs. (13.171) . (13. 8p . 
In each plot, the different lines are for a series of values of pq — Pq, namely pq — pq = 
-10, -5, -1, -0.5, -0.2, -0.1, 0.1, 0.2, 0.5, 2, 5, 7. 



but differ in the non-conformal cases. We study here how the dimensionful prefactors in these 
equations vary with respect to each other when varying jfXq- Let us define the following 
quantity: 



27rmQ /32 
Nf 



(3.29) 



By Taylor expanding for large pip (small ^^Ag), we find the following leading order expres- 
sion for 6: 



6 
5 



1 + 
1 - 



3 3. 3 3. \ 3 iV/ 
8 5 14 y 327r2 

280 J 327r2 N, ^ 



A 



Q 



+ ... 



+ 



< , 

Pq>0 , (3.30) 

where we have defined pq = Pq — Pq. Both for ^Aq ^ or pg — > cxd, one has 5 — > 1, 
recovering the quenched result. This relation, together with the results in figure HJ allows 
one to compute how the tu's would vary with ^Aq if, instead of (13.81) one decided to use 
(EUD, fl3J[5l) as the definition of a;. 



3.6 The holographic a-function in flavored KW models 

As we have anticipated in section [21 the D3-D7 conifold models under study are plagued 
by various UV pathologies. The behavior of the holographic a-function (see eq. (12.41) ). in 
particular, suggests that any UV cutoff we choose to adopt in the analysis of the mesonic 
spectra has to be located below the discontinuity point pa where the function diverges. 
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Figure 5: The holographic a-function in the massless-flavored KW model as a function of 
p - pLp (left) and of (iV^/iVjA = (327rV3)(piP - p)-' (right). 



Fixing the overall constant factor such that a = axw = (27/64) in the unflavored case, 
the function a{p) for the flavored KW models can be taken as: 



a(p) 



27 



where 



H{p) 



~2T 



(3.31) 



(3.32) 



Let us now consider the flavored KW model with massless quarks. The function H{p), 
which in the unflavored case is monotonically increasing as e^'', now has a maximum at 
Pa ~ Plp - 0.270 This behavior strongly affects that of the holographic a-functioii^ as it is 
evident in flgure [H As expected, the holographic a-function has a bad discontinuity at pa, 
where H\pa) = 0. For p < pa, a{p) is positive and increasing with p, whereas for p > pa it 
grows from minus inflnity to zero. 

Far below the Landau pole, the holographic a-function is only slightly varying with 
{Nf/Nc)X and has an almost linear behavior 



a{p) ^ Nl 



27 
64 



+ 2 X l{)-\gl^Nf ) 



(p < Plp) 



(3.33) 



where gprp ~ e'''^''' . This expression can be read as the leading correction to the effective 
number of adjoint plus bifundamental degrees of freedom, due to the internal quark loops 
(see analogous comments for the entropy of the D3-D7 model in flat space in ref. [20]). 

For the non-chiral massive-flavored KW solution, where the sea quarks have a mass which 
is related to p^, relevant plots are shown in flgure [61 The holographic a-function approaches 



Plp- This behavior does not depend 



^H{p) starts from zero at p — > — oo and comes backs to zero at p - 
on the particular choice of integration constant for the warp factor. 

^It is simple to realize that a{p) does not depend on gs and on the integration constant C3. It just depends 
on PLP — p and is proportional to N^. 
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Figure 6: The holographic a-function in the non chiral massive flavored KW modeL We 
have fixed Pq = p^p — 2. On the right, a zoom around Pg — p^p. 

the unfiavored value axw, for p <ti Pq', then it is very slowly increasing with p up to p = pq. 
Slightly above pq the function a(p) starts increasing faster and then it blows up towards a 
point Pa where the previously noted bad discontinuity shows up. In the present case pa is a 
function of pqi it goes to p^p — 0.27, as in the massless fiavored case, for very small sea quark 
masses {plp — Pq ^ oo) and goes to pip for very large masses. For p > Pa the a-function is 
again going from minus infinity to zero. 

In section 15. 2[ we will suggest a possible relation between the behavior of the a-function 
and that of the (^^g functions in the mesonic spectrum. 



4 Meson excitations in the flavored KS models 

This section follows the same structure as the previous one, but we consider solutions in which 
the conifold metric is deformed and there are three-form fiuxes that affect the UV behavior 
of the background. Namely, the quenched set-up is the Klebanov-Strassler solution [2]. We 
discuss mesonic excitations in this background and in its generalizations with unquenched 
massless [1] and massive [8] non-chiral fiavors. The solution with massless chiral fiavors was 
found in [5], but we will not address that case here. 



4.1 Background and excitation equations 

Let us write the metric in the notation of [16]. In Einstein frame it reads (both for the 
quenched and unquenched cases): 

+ a'A^T) (^hj + hl + hj + ^ + -^{h2h (4.1) 

where we have used the definitions in (]A.9|) . There are non-trivial F(3), F(5), if(3) forms, 
which we will not need in the following. We refer the reader to [1] for details. 
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For the fluctuating probe brane, we will consider a supersymmetric embedding first dis- 
cussed in [16], where it was shown to correspond to a massless non-chiral flavor in this set-up. 
Concretely, the embedding is parameterized as 7 = 0, 5 = const, in terms of the variables 
introduced in appendix |Al We discuss in the following mesonic excitations of this probe 
brane. Note that this simple massless embedding cannot be used in the KW cases because it 
would not produce a discrete spectrum, due to the different IR behavior of the backgrounds. 
Oscillations around massive generalizations of this embedding (in the quenched background) 
have been discussed in [21J, but are beyond the scope of the present work. 

Similarly to ( 13. Sp . let us just consider excitations of the gauge field of the form: 

A = a^{T)e'^''i^dx^' + a,(r)e^'=^/ii . (4.2) 
The corresponding second order equations are 

Q = dr {A^ tanh rd^a^) + h a'A^B"^ tanh r , 



,'coth ^ ^ 
= dr{ — r^d^a, 



1 ^ , 1 X 1 tanh ^ 0,0 T 
— dJh-^) + a'B^ coth - 

2 ^ ^ 4 h " 2 



(4.3) 



h 

Again, we want to factor the meson masses as a dimensionful scale related to certain IR 
quantities (evaluated at the tip of the oscillating brane, namely r = 0) times a purely 
numerical factor which we will study. As a first step, let us introduce a similar factorization 
for the function h: 

n-ie*)L„'-^^Hr), (4.4) 

/i3 



and let us define: 



1 



M„,, = 27r2i^ 0^,,, , (4.5) 



where we have introduced the string tension and the 't Hooft couplinj^ in the IR as: 



4 



T, = J-,(e*/'^V^)|„. = ^i|i^. A„ = 8„.MM|„„, (4,6) 

where we have defined Jq = I\t=o- Our goal in the following is to determine the value of 
the w's in the set-ups with and without dynamical flavors. As in the previous section, WKB 
estimates will be very useful. In this case, for both the vector and scalar modes, the integral 
(IB.12P which governs the behavior of the spectrum reads: 

S = ^-3Vy^ 1-2 BdT, C = ^- (4.7) 



^"^In I16j, a constant dilaton was considered but it turns out that a running dilaton does not modify the 
expressions (I4.3p . In the very far UV these equations coincide with those of KW, as expected, see eqs. (13.61) . 
(|3.7p . Looking at the far UV amounts to setting e^^'^'^^P — 0, tanhr — tanh J = coth § 1- One also has 
to identify e^s = GA^, e^-^ = 95^ and r = 3p. 

As in footnote El we emphasize that we abuse of language when calling A the 't Hooft coupling. 
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4.2 Mesons in the quenched approximation 

Let us write the value of the functions that determine the geometry when no backreacting 
flavor branes are present p]. Define: 

Kir) = '-"'"f ' - ^^'^ . (4.8) 
2 3 sinh r 

Then, the different functions of the ansatz read: 



4 



A\T) = ^coshTK{T) , B\r) = -^^, e'^ = const , (4.9) 

whereas the function I is given in terms of the following integral: 

^, , ^ r X cothx — 1 , , , . . 1 , 

/(r) = Jo - / 2 (sinh(2a;) - 2a;)3rfx . (4.10) 

Jo sinh X 

The condition I{oo) = fixes: 

Jo = 0.71805 . (4.11) 
In this case, the second order equations (14. 3 p reduce to: 



= dr {K{t) sinh rS^a^) + < Jq ^/(r 



sinh r 
6ir(r) 



a,; 



„ ,'coth ^ „ 
= a. ( -J^dras ] + 



(4,12) 



The discrete set of u's which make the fluctuations regular in the IR and normalizable in 
the UV can be determined numerically via the shooting technique. We obtain!^ 

co^ ^ 1.367,2.553,3.752,4.956,6.164,... 

tUs ^ 3.100,4.491,5.778,7.032,8.271,... (4.13) 
4.2.1 WKB estimates 

Applying the WKB approximation (as in appendix [B]) to this case, we get the following 
estimates for the masses: 



7 5 

Us ~ C\ln'^ + ^n + -\ogn + 2 , (4.14) 



The same computation has been performed in |16j . Our numerical values differ from those displayed 

_ 1 

there (even taking into account that there is a factor of /q difference between our definition of uj and the 
definition of A in [16j). 
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where ( = 1.214 is what comes from the WKB integral (14. 7p . in which Tcut = oo. The 
constant terms inside the square roots, as well as those with logarithms, do not come from 
WKB but from fitting the points obtained numerically. We have included those terms in 
order to fit masses with low n (the WKB result is only guaranteed to be accurate for large n). 
We stress that such dependences on n are not exact, but just useful simple approximations. 
A comparison of the expressions (I4.14p to the first few numerical values is given in figure [3 




Figure 7: Comparison of the estimates (14.141) to the first few points obtained numerically. 



4.3 Mesons in the massless unquenched set-up 

The deformation of the solution due to a set of smeared massless flavors was found in [A]. In 
the following, we employ a notation slightly different from the one used in [4|. Let us call 
the position of the Landau pole and define: 

^ [2(r, - r)(sinh(2r) - r) + cosh(2r) - 2r r, - 1]^ 
(4rrf) 3 sinh r 

which is such that in the quenched limit lim^^^^^oo A(r) = K{t). The other functions read: 

A'{t) = ^ cosh r A{t) , B\r) = ^^^-^^ , = ^ ^. (4.I6) 

The first two expressions reduce to those in (14. 9 p when — * 00. From the expression of the 
dilaton and the definition of Aq in (14.60 . one gets: 

The remaining function entering the metric is /(r) which can be written as: 

^, , ^ „_i I r a;cothx — 1 — cosh(2x) + 4a;^ — 4x + 1 — (x — 2rd) sinh(2a;) , 

/(r)=/o-2 -AT^ / 2 ^ — ^ ^ — 

Jo [Td - x)^ smh X (cosh(2x) + 2^^ - 4x - 1 - 2{x - r^) sinh(2x))" 



3 



(4.18) 
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which consistently reduces to (14.10p as — oo. We still have to give a condition that fixes 
Jo for a given r^. Notice that, regardless of the value Iq > that we may insert, the function 
/(r) vanishes at some value Th < t^, since, in fact, it is easy to check that I{Td) = — oo. 
However, when ^ 1, there is a region 1 <C r < r^^ — e, where the derivative of /(r) is 

exponentially small. The parameter e ~ rje s^'' is an exponentially small quantity defined 
by the fact that the {x — t^)"^ in the integrand of fl4.18p starts dominating and makes dri 
diverge. We define Iq as the offset of /(r) such that the almost flat region 1 ^ r < — e 
lies around / = 0. For instance, we can set Th = Td — 1 and flx Iq by requiring /(t/j) = 0. 
Any other similar deflnition, as for instance Th = — 2, would only affect Jo by a negligible 
exponentially small quantity (of order e~3'^'"). Thus, as anticipated in section [21 the choice 
of UV cutoffs only affects the IR physical quantities such as meson masses by negligible 
amounts. We emphasize that for this reasoning to be valid, it is necessary that ^ 1, 
making it possible to have a large separation between the UV and IR scales (see related 
comments in section [3.3. ip . 

Let us now write the excitation second order equations (14. 3 p for this case: 



= dAA{T)smhTdra,)+ujll,'^I{T)'^^^!^ a, , (4.19) 

6A(r) Td 



^ , ' coth ^ ^ , 



^^iTi Itanhf 1 Td-T T 



The values of lu only depend on which, in turn, is related through (I4.17P to ^Aq, the 
physical parameter controlling the deformation produced by the unquenched flavors. 

Notice that the D3-D7 background in [1] has a (good) curvature singularity at r = 0, the 
minimal distance reached by the D7-brane probe we are focusing on. Nevertheless, since 
the metric components are not singular there, we can extrapolate the fluctuations up to the 
origin where, as usual, we impose them to be regular. 



4.3.1 WKB estimates 

Similarly to section [3l the WKB formulae (I4.14p also match the results in unquenched back- 
grounds quite well, the main effect of the flavor backreaction being a shift in the value of (. 
Inserting fl4.16p in (14. 7p yields: 

\io v^A(t) 

We will choose Tcut < ^ Th — 0.38. It is worth stressing once again the independence of 
the result in such UV prescription: for 1 <^ r < r/^, the integrand is of order e~'^^^ since 
/ oc e~3'^ and A oc e^^^^. Had we chosen a different Tcut, it would only have changed the 
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integral by an exponentially small quantity. In more physical terms, the main contribution to 
the integral which fixes the meson masses comes from the IR region, as expected on general 
grounds. Following the reasoning at the end of section I3.3.2[ one can check that the size of 
the corrections associated to the UV features is of order 4^- 

J^vv 

By expanding for large r^, we find: 

C ^ 1.214(1 - 1.14 X lO^^^Ao + ...). (4.21) 

4.4 Mesons in the presence of massive dynamical quarks 

We now consider mesons in the fiavored KS solution with massive dynamical quarks [S]. 

Let us begin from the step function approximation for the function representing the ef- 
fective number of fiavor degrees of freedom, A^/(r) = A^/6[r — r,], where Tq is related to the 
mass of the dynamical quarks. The solution in the IR, for r < r^, is just the ordinary KS 
solution but for the value of Jq, which is again set by the requirement that hijh) = 0. For 
T > Tq the solution is practically equal to the massless fiavored one fJj considered in the 
previous section. The only difference is in the integration constant, previously set to one, in 
the function A: 



[2(r, - r)(sinh(2r) - r) + cosh(2r) - 2r r, - (cosh(2r,) - 2r^)^ 

= ,1 ■ , > (4-22) 

{4:Td) 3 smh T 



and in a relation between the IR and UV scales /i: 

4/3 _ 4/3 f Td 



1/3 



These relations follow from the requirement of continuity of the metric at Tq. Also, the 
coupling in the IR will now read: 

Ao = 87rg,Me%=o = A., = Jf""'^ . . (4.24) 

Nf{Td-rq) 

It does not run below the sea quark mass scale Tq. Using the definition (14. 4p we have, in the 
UV region]^ 



Hr) = -2-hliT,-rq)J 



X cothx — 1 



Td — sinh^ X 
— cosh(2x) + 4x^ — 4x + 1 — (x — 2Td) sinh(2x) 
(cosh(2a;) + 2x'^ -Axra- (cosh(2rg) - 2t^) - 2{x - r^) sinh(2x))^ 



dx . (4.25) 



""^^Note that the continuity of the metric at Tq impHes a discontinuity of /: luvijq) = IiF/.iTq)(^ 



2/3 



Td-Tg 
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The equations for the mesons are the same as in the previous sections and their mass 
is controlled in the WKB approximation by the function S. It is calculated with the KS 
solution for t < Tq, while for r > it reads: 



-•uv 



The result for the function 



Tcut T 4 T 

-'o ^ 



V) 2 



V6A(r) 



1/6 



dr . 



c = - 



(4.26) 



(4.27) 



which gives the meson mass dependence on the number of flavors and their mass, is shown 
in figure [HI The behavior of the masses is qualitatively the same as for the KW model, 




M 



Figure 8: The plot of (, see (14.271) . entering the WKB meson mass formulae. The different 
lines are, from bottom to top, for Tg = 0.1, 1, 5, 10, 15, 20. The meson masses vary only very 



Nf 



mildly when changing 
large quark mass. 



An 



Also, they go to the quenched case both for small jjXo and for 



figure m The masses are mildly, almost linearly decreasing with -jjXq. Moreover, the more 
massive are the dynamical quarks, the smaller is the variation; for very large quark masses, 
the meson masses are practically constant and equal to the quenched values, as expected. 

We have repeated this calculation for the exact solution found in [8], where the function 
Nf{T) is not of the approximate step form but is explicitly calculated from the brane em- 
beddings. Since the formulae for this background are rather involved, we do not report them 
here and refer the reader to [8j. In any case, the qualitative behavior of the calculation of 
the meson masses is exactly as described above in the step approximation for Nf{T). The 
plot of the function ( is qualitatively identical to figure [H also in the region not accessible 
by the step approximation, where the flavor branes reach the tip of the conifold and so the 
sea quark mass is smaller than the dynamically generated scale, but it is still non-zero. 
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Figure 9: The holographic a-function in the massless-flavored KS model. We have set 
Td = 10,/i(0) = 1 and fixed 2^/^gsi32TT^)M'^ = Snii^'^N f^. In this case Th ^ 3.93 and 
r, ^ 3.57. 



4.5 The holographic a-function 

The expression for the holographic a-function in the fiavored KS case is given by: 

(r, - r) h{r) 



(3{t) = a' 



6 r, A(r)2 ' 



3 Td 

a(r) ~ 21(3'^'^H'"^ [H'] . (4.28) 

It is clear that the holographic a-charge can be real only for r < t^, where h{Th) = 0. For 
r — > the function goes to zero as 

a(r) ^ Xl II r^ (r ^ 0) , (4.29) 

and increases with r up to the point < Th (where H'{Ta) = 0) where it diverges and it 
has the bad discontinuity discussed in section [21 A representative plot of the holographic 
function a(r) in the massless-fiavored KS solution is given in figure [H 



4.6 The domain wall tension 



We compute here an expression for the domain wall tension, which is related to the value of 
the gluino condensate. It comes from the action of a D5-brane sitting at r = 0, wrapping 
the finite [22], namely: 



1 



T=o gs 27r'^a'^ 
We can write it in terms of the IR quantities defined in (14.6 

Tdw 1 



' 62 



87r2v^ 



3 
-In 



r=0 



(4.30) 



(4.31) 
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This expression gives a meaning to Iq in terms of physical quantities. By expanding for large 
Td, one can find the effect of the (massless) unquenched quarks on this quantity to be: 

Iq ^ 0.71805(1 + 1.16 X 10-2-1 Ao + ...). (4.32) 

5 Towards an interpretation of the behavior of the 
mass spectra 

In the previous sections we have extracted the mass spectra for specific mesons in the KW 
and KS models with and without dynamical sea quarks, by studying fluctuations of the 
U{1) gauge field on the worldvolume of a probe D7-brane. The equations of motion for more 
generic fluctuations can be very difficult to solve and/or diagonalize and so it could be useful 
to have an idea of how the mass spectra behave (as a function of the scales and couplings 
of the theory) in those cases. Moreover it would be nice to understand the behaviors we 
have found for the mesonic spectra at fixed scales or couplings in terms of simple general 
expressions. 



5.1 Unflavored KW-like models 

Let us start by reviewing some known results in the unflavored D3-brane models on the 
singular conifold, or, in more generality, on a singular toric Calabi-Yau cone over a 5d 
Sasaki- Einstein compact manifold X5. The gauge theory on the D3-branes has gauge group 
of the form SU{Ny where p is an integer, and there are bifundamental and adjoint matter 
fields coupled through various superpotential terms. The dual supergravity background has 
a near horizon AdS^ x X5 metric: 

ds^ = —dxjx^' + — idr"^ + r'^dsVj , (5.1) 

with AdS radius given by: 

This background is dual to an IR surface of fixed points for the planar gauge theory at strong 
(effective) 't Hooft coupling A = QprpN^ ~ QsNc ^ 1, where gi = g-z = ■■■ = gp = gpT- The 
volume of X5 is holographically related to the IR central charge apT as in 



Let us now place a (supersymmetric) probe D7-brane on AdS^ x X5, assuming that the brane 
reaches a minimal distance tq from the origin of the space transverse to the D3-branes. The 
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D3-D7 open strings will provide fundamental flavor flelds of mass tuq = tq/ {2Tia'). The low 
spin (J = 0, 1) mesons dual to fluctuating worldvolume flelds on the D7-branes will have, 
just as in the simplest X5 = case |2l], masses of the order: 




M~^~^,/^, (6.4) 

where we omit universal (i.e. independent on X5) numerical factors. Using the relation 
27ia'mQ = tq = R^y Gtt{rQ) we can easily rewrite the formula above as: 

, (5.5) 



AV4 \aFT, 

where 2Tra'TQ = Guijo). This expression reproduces (in the constant Xq = \ case) the 
parameterization of the meson masses introduced in eq. ( 13. Sp . once adapted to the conifold 
case, where X5 = T^'^ with Vol{T^'^) = IQ-k^/27 and so apx = N^{27/64:). 

As already outlined in the literature, these results imply that at strong coupling the 
binding energy of the QQ system is very large, since the mesonic masses are parametrically 
smaller that the total quark mass 2mQ. Moreover, as was observed in [25] for the AdS:^ x 
case, they imply that the typical size of the mesons is L ~ ^/X/niQ which is very different 
from the weak coupling expression (the Bohr radius) L^^ak ~ ("^q-^)"^ and much larger than 

For our purposes, a formula like fl5.5p allows us to show how the meson masses are ex- 
pected to scale with the effective number of (adjoint+bifundamental) degrees of freedom 
(which are accounted for by apx)- For flxed Tq/^/X the masses decrease as this number is 
increased. This behavior can be understood as an effect of the increasing binding between 
the constituent quarks as apT is increased. 



5.2 Flavored KW models 

It is natural to think about a possible extension of formula ( 15. 5p to the case where the probe 
D7-brane is put in the fully backreacted D3-D7 background dual to, say, the flavored KW 
model. Since the fluctuating flelds which are dual to the low spin mesons are localized on 
the probe D7-brane at pq (we use, as in the rest of the paper, a different radial variable 
when referring to the unquenched solutions) we can expect the mesonic masses to scale as: 



^ (5-6) 



where Aq ~ 'N^.e'^^^'^^ and aq = a^pg) are the 't Hooft coupling and the holographic a-function 
(see eq. f l3.3ip ) evaluated at pg. The expression above could provide a simple interpretation 
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Figure 10: The function (q f lS.Sp in the massless flavored KW model. 



of the linearly varying (with NjXq/Nc) functions (s,v entering the mass parameterization 
(ESD: 

Cs,JsAn)- (5.7) 



If we define: 

Co=( 

we can first notice that, just as (s,v, this function is indeed always decreasing (outside the 
pathological UV region) as NjXq/Nc increases. In fact, as we have observed in section [3^ 
the holographic a-function is always increasing slightly with g^j^Nj ~ NfX/Nc as a result 
of the increasing number of effective (adjoint plus bifundamental) degrees of freedom due to 
the internal quark loops. 

Moreover, in the massless-flavored KW case, we can see that the function (q also has a 
quantitative behavior which is close to that of (s,v Far below the Landau pole, (q is in fact 
a linearly decreasing function of {Nf/Nc)XQ (see figure [TOj) . The linear behavior persists up 
to large values (0(200)) of {Nj /Nc)Xq. The best linear fit is actually: 

where the term in parentheses is quite close to what we have found for (s,v in eqns. fl3.22p . 

Analogous comments apply to the massive-flavored KW model, though the quantitative 
matching between (q and (s,v is realized only in the pg > Pg region. 

A proposal for an alternative way of interpreting the behavior is briefly discussed in 
appendix O 



5.3 Flavored KS models 

A parameterization of the meson masses analogous to the one in (15. 7p has been adopted in 
section m in the (flavored) KS cases. The only differences are that: 1) in those cases we have 
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just focused on mesons made up of massless constituent quarks (the corresponding probe 
D7-brane reaching the origin); 2) the effective 't Hooft couphng was written in terms of the 
number M of fractional D3-branes instead of that {Nc) counting the regular D3-branes. The 
first point just implies replacing Tq with Tq = Tg and the second just means that we have to 
replace ^/\Q with Aq in the mass formula. The resulting parameterization is that written in 
eq. (14. 5p . In order to give an interpretation to the C, functions in the flavored KS cases one 
could thus think of naively extrapolating the previous results and so to relate those functions 
to , with ao being the holographic a-charge at the origin. However this reasoning would 
certainly fail, since = in the KS models. 

For the cases where mg is much larger than the IR dynamical scale, we expect the qual- 
itative considerations proposed for the flavored KW models to apply. Instead, we should 
consider alternative ways to interpret the C, functions in the mq = flavored KS cases. One 
interpretation is the following. From the KS-like metrics we read that Ts/m\^ ~ XqI^^ , 
where ttikk is the mass scale of the glueball and KK modes. Since in our analysis of the 
spectrum we keep fixed the ratio (T^/ Aq)^''^, we are keeping fixed iriKKiy^ ■, which means that 
fnKK ~ Iq^^^- If "we compare formulas f l4.2ip and f l4.32p we are led to argue that C, ~ Iq^^^ is 
a measure of the variation of the KK (or glueball) mass scale with the number of dynamical 
sea flavors and their masses. Possibly this variation can also be related to a change in the 
effective number of degrees of freedom around r = 0. From eq. (14.291) . in fact, we get that 
[a(r)/(M^r^)]~^/'^ ~ X^^^'^Iq^^'^ which displays a qualitative behavior similar to that related 
with rriKK- 

6 Summary and outlook 

Since in real world QCD the number of flavors is of the same order of the number of colors, it 
is of obvious interest to study the extensions of the standard string/gauge theory correspon- 
dence to theories including the dynamical effects of fields in the fundamental representation 
[HI EH [12]. Some technical difficulties in finding such extensions can be overcome if the 
branes carrying flavor degrees of freedom are smeared in the internal directions |271 [13] . De- 
spite the fact that, due to this simplification, there are nowadays a large number of available 
models in the literature [3]- [8], [T3l [131 [T5] . some of their relevant features have not been fully 
investigated. Among the others, an important issue concerns holography in the presence of 
a UV cutoff. In this paper, we have addressed some relevant topics in the flavored versions 
of the Klebanov-Witten [3], [7] and Klebanov-Strassler theories [H [8], dual to D3-D7 brane 
systems on the (deformed) conifold. 

We have started by refining the range of validity of the supergravity solutions. Besides the 
point where the dilaton diverges, which is commonly associated to the presence of a Landau 
pole in the dual field theory, we uncovered in these backgrounds a point of singularity in the 
holographic a-function. The existence of the singular point, which is located in the UV of 
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the solution at least for small enough Nf, is clearly related to the presence of D7 branes. 
The presence of the Landau pole already required a cutoff. The new effect described here 
shifts the needed cutoff towards smaller energies. 

Keeping in mind this issue, we have performed a study of small spin mesons in our conifold 
models. The meson masses, even in the quenched case, are quite different from the QCD 
ones, due to the fact that the mesons are very tightly bound [24]. So, there is no direct 
way of matching the corrections due to unquenching with actual QCD data. Nevertheless, 
our perspective is that the studies in the calculable gravity duals can be useful to gain some 
qualitative insights on the way (and to which extent) the screening effects can be understood 
from simple considerations on the physics of the models at hand. 

Concretely, we have introduced a probe D7-brane in the backgrounds and have calculated 
the masses of some vector world-volume modes. Our primary interest was to investigate 
how the meson masses depend on the dynamical flavor degrees of freedom, in particular 
their number and masses. The presence of fields in the fundamental representation screens 
the charges, so the meson binding energy is affected by flavors. The comparison of the 
screening effect among theories with different number and mass of the flavors requires some 
conventional choice of the physical scale to be kept fixed in the process. Once this choice 
has been performed, the spectra of mesons have been calculated numerically, both for the 
quenched and the unquenched KW and KS theories. Keeping fixed the coupling at the probe 
quark scale yields a decreasing of the meson masses as g'px^f is increased. Had we kept 
fixed the coupling at some larger UV scale, larger g^-pNf could have yielded larger meson 
masses. 

Our computation of the change in the meson masses due to the unquenched flavors is, in 
spirit, equivalent to the Lamb shift in atomic spectra: both are effects that show up when 
taking into account the (slow) running of the couplings in a theory with a Landau pole. The 
main qualitative difference is that while in QED the theory is weakly coupled in the IR, in 
our case the theory is always strongly coupled. 

Finally, we have suggested that the mass dependence on the flavor parameters in the KW 
models seems to be simply described by the behavior of the holographic a-function, counting 
the number of effective degrees of freedom at a certain scale. In the KS cases it is more 
plausibly related to the changes in the glueball (or KK) mass scale as the dynamical flavor 

^''A possible explanation for such pathology is the breakdown of the naive interpretation of the gravity/field 
theory duality, due to the incomplete freezing of the gauge coupling on the D7 branes. In fact, if the 
D7 SU{Nf) (or, better, U{l)^f) is not really a global symmetry but rather a (very) weakly gauged one, 
at a certain UV energy scale its gauge dynamics will be dominant over the D3 one. At that point, the 
interpretation of the system as a 4d gauge theory coupled to Nf flavors breaks down. Thus, the latter 
interpretation is trustworthy only up to some UV cutoff scale. An alternative explanation is that in the UV 
the claimed decoupling of the gauge theory with the bulk gravity modes is not fully realized. The presence 
of the singularity in the UV prevents us from casting any reasonable conclusion. Wc thank Carlos Nunez 
and Francesco Benini for related comments. 
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parameters are varied. 

A relevant case that remained out of our investigation is the spectrum of mesons for massive 
probe quarks in the (flavored) KS model. In general, the study of other observables in such 
flavored backgrounds could shed further light on the holographic correspondence in these 
models with UV cutoffs. It would be interesting, in particular, to investigate the behavior 
of the entanglement entropy. Finally, it would be important to understand holographic 
renormalization in theories with UV cutoffs. 

Acknowledgments 

We are grateful to R. Argurio, F. Benini, C. Nunez and J. Shock for relevant comments and 
to E. Imeroni and S. Kuperstein for useful discussions. This work has been supported by 
the European Commission FP6 programme MRTN-CT-2004 v-005104, "Constituents, fun- 
damental forces and symmetries in the universe". F. B. is also supported by the Belgian 
Fonds de la Recherche Fondamentale Collective (grant 2.4655.07), by the Belgian Institut In- 
teruniversitaire des Sciences Nucleaires (grant 4.4505.86) and the Interuniversity Attraction 
Poles Programme (Belgian Science Policy). A. C. is also supported by the FWO - Vlaan- 
deren, project G. 0235. 05 and by the Federal Office for Scientific, Technical and Cultural 
Affairs through the Interuniversity Attraction Poles Programme (Belgian Science Policy) 
P6/11-P. A. P. is also supported by a NWO VIDI grant 016.069.313 and by INTAS contract 
03-51-6346. A. R. is also supported by the MCINN and FEDER (grant FPA2008-01838), the 
Spanish Consohder-Ingenio 2010 Programme CPAN (CSD2007-00042) and Xunta de Galicia 
(Conselleria de Educacion and grant PGIDIT06PXIB206185PR). 

F. B. and A. L. C. would like to thank the Italian students, parents and scientists for their 
activity in support of public education and research. 



A Some technical details of section 3 



A.l The probe brane embedding 

In order to study mesonic excitations, we consider a flavor D7-brane which preserves su- 
persymmetry in the backgrounds (13.11) . In particular, we take a probe brane embedding 
Z4 — Z3 = fi which introduces non-chiral flavors and was flrst studied in detail by Kuperstein 
in [T6] . In the following, we follow the steps of [16] in order to introduce a new set of variables 
that will be convenient to describe this embedding and its excitations. The (dimensionless) 
z's are defined as: 

= el PgtC'^-^i-^^) sin ^ sin ^ , Z2 = e^" e^^^^^^'^^'^ cos - cos - , 

^ 2 2 2 2 

^3 = ei''e^(^+^^-^2) cos - sin - , ^4 = e5''e5('^-^i+^2) sin - cos - . (A.l) 

2 2 2 2 
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The conifold equation ZxZ2 = can be written as: 

detW = , with W = ( ] . (A.2) 

\ —Zi Z4 

In order to write it in terms of the angles, let us introduce SU (2) matrices: 

_, cos 1 63 ('^'+^') -sin|e-5('^-^>) 

sin|e5(V'«-^») cos|e-5('^»+^0 ' ' ^^'^^ 

Since the ipi only appear in the combination ipi + let us set V'l = "^2 = ^i^- We can 
explicitly check from flA.ip . flA.2p that: 

W = el'>L,(^l (A.4) 

Now, the trick of [16] consists of introducing a change of coordinates which eliminates 62, v?2 
in favor of two new angles 7, S with ranges 7 G [0, vr], 5 G [0, 47r), as: 

L,S, with 5=1 "''^^i^e;.^' 1 . (A.5) 

-ism^e a cos 2 

The nice point is that this new system of coordinates is very well adapted to the massless 
non-chiral embedding TiW = Z4 — z^ = fJ^, since by explicit computation we find that: 

TiW = i sin ^e'h^f = /i . (A.6) 

3 

Suppose without loss of generality that /i is real and positive so we can define /i = e^PQ. 
Then, the Z4 — z^ = fi embedding is just: 

sin |et^ = et^Q , (5 = Svr . (A.7) 

This embedding equation is valid both for the quenched and unquenched set-ups. In par- 
ticular, in the unquenched case, it is possible to show that flA.7p gives a solution for the 
worldvolume action Sdy = T-j j {—e'^ — <lei{P[G + -P[C'(8)]), where the P denotes the 
pull-back on the brane worldvolume. To do this, it is necessary to use (13.41) and take into 
account that C(8) (defined as dC{^) = e^'^(*F(9))) takes a non-trivial value in the background: 

QsNf 6^*^+^^ 

C(8) = — — sin 61 cos 7 {dxo A dxi A dx2 A dx^ A dp A dOi A d(pi A dip) . 

(A.8) 
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A. 2 The meson excitation equations 

In order to study the fluctuations, we want to write down the metric in terms of the new 
coordinates 7, 6. Let us start by rewriting the original metric in terms of two sets of left 
invariant one forms: 

hi = — cos sin 9id(pi + sin ^dOi , h2 = — sin sin Oidipi — cos ^dOi , 

hi = — cos sin 92dip2 + sin ^^62 , h2 = — sin ^ sin 92dip2 — cos ^^62 , 

/is = ^ + cos^ifiv^i , /is = + cos92d(p2 ■ (A. 9) 

The angular part of the metric in (13 .ip reads: 

d4 = ^(hl + hl + hl + hi) + ^(/is + /is)' . (A.IO) 

Substituting the expression above we get the usual metric written in terms of the angles 
Oi) fi, ^2, V^2, The form of the /ij's in terms of hi, 7, 6 can be obtained by explicit compu- 
tation from (1X51) [16]: 

/ii = (/ii — (^7) cos5 — (/is sin7 + /i2 COS7) sin 5 , 
/i2 = (/ii — (^7) sin5 + (/i3 sin7 + /i2 C0S7) cos5 , 
/is = (/is cos 7 — /i2 sin 7) + . (A. 11) 

Substituting in flAlOl) we get: 



dsl = ^ I ^^7^ + 2 f /ii - Y ) +/i^ + (sin7/i3 + cos7/i2)' I + 



2 



e 



2/ 



+ ^(d5 + /is(l + cos7) -/i2 sin 7)2 . (A.12) 
9 



The induced worldvolume metric is obtained by inserting flA.7p in this expression. In order to 



compute the meson spectrum for the modes discussed, we insert the ansatz for the fluctuation 
(13. 5p in the equations of motion stemming from the worldvolume action: 



Sd7 = T^[- j d^xe^^Jdei{-P\G^,]+2^la'e--2F^,) 

+ / P[Ci.)] + / ^[^4)] A F A f) . (A.13) 

Keeping only linear terms in the fluctuation equation, one flnds (13. 6p . (13. 7p . 
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B Schrodinger form and WKB estimates 



In this appendix we shall start by reviewing the mapping of a generic mesonic fluctuation 
equation to a Schrodinger wave equation and the subsequent use of the WKB approximation 
to get an estimate of the corresponding mass levels. We shall follow closely section 3 of ref. 
[28] , which we will adapt to our notations. Let us suppose that the fluctuation is represented 
by a function 0(p), which satisfies a differential equation of the form: 

e-f^d^{e-f^F{p)dp4>) + {^'H{p) + P{p)) 4> = , (B.l) 

where the radial variable p is non-negative (p > 0), is a mass parameter and F{p), H{p) 
and P{p) are three arbitrary functions that are independent of uj. We will assume that p 
takes values in the range < p < pcut, where pcut is large and, in the quenched cases, we 
will just take pcut = +oo. 

By performing a suitable change of variables, eq. (IB. 11) can be written as the zero energy 
Schrodinger equation: 

d^i: - V{y)ij = 0, (B.2) 

where V{y) is some potential to be determined. The change of variables needed to pass from 
f lRTl) to (El is [28]: 

y 

6^ = 6^-1, (t) = -^ip . (B.3) 



Notice that y > — oo. In order to write the expression for the potential V{y), let us define: 

Fo = e-yp , Ho = e^H , Pq = P . (B.4) 
Then, the Schrodinger potential V becomes [28] : 



where the primes denote derivatives with respect to y. Notice that the potential V in fIB.SP 
depends parametrically on uj. When V satisfies some general conditions, uj can be fine-tuned 
to some particular discrete set of values in order to produce a normalizable zero-energy 
solution of the Schrodinger equation (]B.2p . 



One can get rather accurate estimates of the mass levels by means of the WKB approxima- 
tion, whose starting point is the standard semiclassical quantization rule for the Schrodinger 
equation, namely: 

{n--)7T=l dy ^/-Viv) , n>l, (B.6) 
yi 



2' 



where n E 21 and yi and y2 are the turning points of the potential (y{yi) = V{y2) = 0). One 
can evaluate the right-hand side of eq. ( IB. 61) by expanding it as a power series in l/u. By 
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keeping the leading and subleading terms of this expansion, one can obtain as a function 
of the quantum number n. Actually, as shown in [28], one can obtain the WKB levels by 
simply looking at the IR (p — > 0) and UV (p oo) behavior of the functions F, H and 
P entering the original fluctuation equation ( IB.ip . Let us assume that near p 0, oo these 
functions behave as: 



Fip'' , H ^ Hip'^ , P ^ Pip'' , as p ^ , 



F^F2 e^'P , H e^^P , P ^ P2 e''^^ , as p ^ 00 , (B.7) 

where Fj, ifj, Pj, Si and are constants. The consistency of the WKB approximation 
requires [28] that S2 — Si + 2 and ri — r2 — 2 are strictly positive numbers, whereas S3 — Si + 2 
and ri — rs — 2 can be either positive or zero. Let us define, following ref . [28] , the quantities 

a\ and /3i as: 

ai = S2 - si + 2 , A = ri - r2 - 2 , (B.8) 

and 0:2 as: 

ifs3-Si + 2^0, 

a2 = { (B.9) 

si - 1)2 - ^ if S3 _ + 2 = . 

Similarly, we define ^2 in the form: 

|ri — 1 1 if ri — r3 — 2 7^ , 

/52 = { (B.IO) 

y/(ri - 1)2 - § ifri-r3-2 = 0. 

Notice that a\^2 are determined by the IR behavior of the fluctuation equation, whereas /3i^2 
can be extracted from the UV limit of F, H and G. Actually, the mass levels for large 
quantum number n can be written in terms of a\^2 and /?i^2 as [28j : 




^WKB = ^ ^ - 1 + — + ^ h > 1) , (B.ll 



where S is the following integral: 



pent 







*e^/||. (B.12) 
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B.l WKB for KW backgrounds 



Let us now particularize our formalism to the general case of a metric parameterized by two 
functions / and g and a warp factor h as the one written in (13.11) . The differential equations 
for the vector and scalar modes (13. 5p have been written in eqs. (13. 6p and (13. 7p respectively. 
Let us rewrite these equations in terms of the variable p = p — Pq defined in (13.110 . One 
has: 



l + e 



1) a, = , (B.13) 



4 ie^P -1)1 



a., = 



:B.14) 



where /, g should be considered as functions of p, u^^s have been defined in (13. Sp and / and 
Iq are given by: 



lip) 



h{p) 



Hp = Pq) ■ 



(B.15) 



277rgsN, 

By comparing the fluctuation equation ( ]B.13p with our general formula ( IB.ip . it follows that 
the integral determining the mass gap for the vector modes is: 



Pcut 



,2/(p)+3p _^ |g29{p) _ g2/(p) j^-^^ 



d3p 



:B.16) 



Similarly, for the scalar modes in f lB.14p the corresponding S-integral is: 

Pcut p2/(p)+3p _^ I g2g(p) _ g2/(p) j^-^ 



d3p 



(B.17) 



B.1.1 Vector modes in the quenched KW 

Let us apply our previous formalism to the case in which the background is just the KW 
without including the backreaction of the flavor branes. In this case we must simply take 
/ = fi' = P + Pq s-iid / = e~^^~^'''3 in our general formulas (see (I3.10p ). One can check that, 



in this case, the fluctuation equations flB.13P and flB.MD reduce to the ones written in (13.120 
and (I3.13p . Notice that in this case the range of p is just < p < +oo. After removing a 



common constant factor the functions F, H and P take the form: 



o}-p _ 1 



-3p 



1 - 



P{p) = . 



(B.18) 
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Changing variables as in (1B.3P we can get the expression of the Schrodinger potential (IB.Sp 
for this case. One gets an expression of the type: 



A{ey + iY{e^y + ?>ey + 3)2 

where V^iy) is a polynomial in of degree seven, namely: 

V^{y) = Ae'^y + 34e^y + {129 - Aul) e^^ + (274 - 24cj^) e^^ + {3A6 - 60ujl)e^y + 



B.19) 



+ (258 - 75ujl)e^y + 15(7 - 300^6^ + 9 (2 



It follows from the above expression that: 



lim VM = 



lim K(l/) 



1 . 



y^oo 



(B.20) 



(B.2i: 



In general, when is above some threshold (related to the meson mass gap), the potential 
V has a minimum around y = 0, where is negative (see figure [TTj) . 





Figure 11: Schrodinger potentials for the vector (left) and scalar (right) excitations in the 
quenched KW model. We used = 3 and = 6- 

Let us now estimate the WKB energy levels for the quenched vector modes. By analyzing 
the behavior of the functions flB.lSP near p ^ 0, we get that the coefficients defined in flB.7p 
are: 

3 

Fi = 3 , si = 1 , Hi = - , S2 = , 



while Pi = S3 = 0. Moreover, the behavior at p 00 is characterized by: 

F2 = 1, n = 3, H2 = I, r2 = -1 , 

with P2 = = 0. Using these values in eqs. fIB.SP - flB.lOp we get: 



(B.22) 



;B.23) 



02 = , 



A = 2, 



/?2 = 2, 



;B.24) 
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and then, as a2/«i + /?2/A = 1, the WKB formula for the masses becomes: 

TT 

uj^ = — n , (B.25) 

with: 



i r ,,,-.v^ ^ f$(lF( - i, i. i. i) . (B.26) 



Numerically ^ 1.3285, which, after substituting this value in flB.25p . yields the expression 
of displayed in (13.171) . 

B.1.2 Scalar modes in the quenched KW 



Similarly, by using the functions /, g and h of (13.101) in (]B.14p . we can study the scalar 
fluctuations in the quenched KW case. We get: 



F{p) = e'-^ (e^^ - 1) , H{p) = ^^^±i , P{p) = -Qe^'^ 



3 p3p 



(B.27) 

where, again, we have removed an irrelevant constant factor. The potential Vs associated to 
these functions takes the form: 

= 4(e»+ine-' + 3e. + 3) ' ''''■^^^ 

where now the polynomial Vs{y) is: 

Vs{y) = 49e^y + 247 e^^ + {501 - 4uj^^) e^^ + (514 - 12cj,2) e^^ + 

+ (271 - 12ujI) + (63 - Qoo^^) + 3 . (B.29) 

Then: 

1 49 

hm Vs{y) = - , lim = -r ■ (B.30) 

J/— >— oo 4 3/— ►oo 4 

In figure [TT] we have plotted Vg for some typical value of ujs- As before, when Ug is above 
the mass gap, Vg has a unique minimum, where the potential is negative. 

In order the get the WKB estimate of the mass levels, let us look at the behavior of the 
functions F, H and P at p = 0. We get: 

Fi = 3 , si = 1 , ^1 = ^ ' S2 = 0, ^1 = ' S3 = -1 , (B.31) 
whereas, by looking at p ^ oo one concludes that: 

33 

F2 = 1 , n = 5 , H2 = 1, r2 = 1 , P2 = -— , rs = 3 . (B.32) 
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Using these values, one arrives at the following coefficients: 

ai = 1 , a2 = 1 , A = 2 , (32 = 7 . (B.33) 
Notice that P2 must be computed from the square root in (IB.lOp . which gives the integer 



value 7. Moreover, in this case: 

' 4 



1 r,-^.,VTO_ ._1_1_5__1, (^^3^^ 



V2 Jo ^fW^ r(^|) ^ 2' 3' 6' 2 

Numerically ^ 1.53173 and the WKB mass formula becomes: 



a;, = 2.051 \\n[n^^-^ . (B.35) 



2, 

To improve the agreement of this WKB estimate with the numerical results for small n we 
can add a new term independent of n inside the square root in such a way that the numerical 
values are fitted. It turns out that an excellent fit is obtained when this new term is 15/8. 
The resulting expression of lOg is just the one written in (13.171) . 



C An alternative way to interpret the (^'s 

Let us comment on a possible alternative way to interpret the C functions entering the 
expression for the mesonic masses. Let us write the mesonic mass formula in the following 
way: 

M~V^, (C.l) 
Ai(/i) 

where, taking inspiration from the Lamb shift calculations in QED, we take the energy 
scale /i to be related to the (modulus of the) binding energy E\j of the meson, defined by 
M = 2mQ — El,. Let us take = \Ei,\/2 for example. Since M = e2mQ, with e <^ 1 related 
to the inverse 't Hooft coupling (at the scale mq), we have: 

Eb^2mQ{l-e). (C.2) 

This relation suggests to relate the energy scale fi with a radial variable: 

p^PQ(l-eQ), (C.3) 

and so to get, by a simple expansion of (IC.ip . a formal relation like: 

Mr^y^ll + ewipQ)] , (C.4) 

where the function w{pq) has to be determined case by case. The above term in parentheses 
could possibly provide an alternative interpretation of the ( functions. We leave this issue 
for future studies. 
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D Comments on the quark- ant iquark potential and high 
spin mesons 

Let us consider the holographic expression for the potential V{L) felt by a QQ pair probing 
our D3-D7 models. Obviously, due to the presence of UV cutoffs in our set-ups, we cannot 
take these sources as being strictly static: their mass mq, in any case, will have to be smaller 
than the UV scale. 

The QQ meson is dual to an open string attached to the probe D7-brane (at pg), bending 
in the bulk up to a minimal radial position pup. The Minkowski separation L between the 
test quarks, as well as the potential V{L) (i.e. the total energy renormalized by subtraction 
of the static quark masses) depend on pup. For an open string embedding given hj t = T,y = 
(T,p = piy) where y G [— L/2,L/2] is one of the spatial Minkowski directions, one finds [29]: 

GPtip 



:D.r 



L{pup) = 2 / , dp 



V{pup) - ^ 



27ra' 



PUP P^P^ - P4 

PQ Qp rPQ 

dp — I G dp 



PupJP^-P^l J-oo 



where P, G are expressed in terms of the string frame metric (we do not consider BH back- 
grounds so Gtt = Gyy): 

P = Gtt, C = a/ GttGpp , (D.2) 

and the "tip" subindex means that the quantity is evaluated at p = pup. From the relations 
above it is clear that the short distance behavior of V{L) is not Coulomb-like. Instead, at 
L = 0, i.e. at pup = Pq, we get V{0) = -2mQ, where: 

Expanding the previous expressions around L = one actually gets: 

(2,771 \ 

i.e. a linear potential, where the effective string tension Tq is given by 2TTa'TQ = Guipq)- 
This result is generic and only depends on the fact that we are considering a D7-brane at 
finite pq: indeed it was previously noticed for a probe brane on AdS^ x in |2^ . 

The occurrence of this linear potential at small L explains why in [24] the masses of 
mesons with spin 1 <^ J <^ \/~\ follow linear Regge trajectories. The very same behavior 
is thus expected for the analogous mesons in the (flavored) KW models. The Regge slope 
dependence on the sea flavor parameters will be simply accounted for by the expression for 
Tq. The study of the screening effects on this slope will thus easily follow. 
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Mesons with spin J ^ should behave instead as non-relativistic quark- ant iquark pairs 
(also at strong coupling) interacting through the large L limit of V{L) which, in the KW- 
like set-ups, is Coulomb-like. The mesonic masses in the flavored KW cases in such regime 
should have expressions analogous to those in |24j and the study of screening effects on these 
mesons should not be a difficult task. 
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